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The energy loss of heavy quarks in a quark-gluon plasma of finite size is studied within the light-cone 
path integral approach. A simple analytical formulation of the radiative energy loss of heavy quarks 
is derived. This provides a convenient way to quantitatively estimate the quark energy loss. Our 
results show that if the energy of a heavy quark is much larger than its mass, the radiative energy 
loss approaches the radiative energy loss of light quarks. 
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In high energy heavy ion collisions hard scattering of 
partons occur in the early stages of the reaction, well 
before a quark-gluon plasma(QGP) might have been 
formed, producing fast partons that propagate through 
the hot and dense medium and lose their energy. Hard 
hadronic probes have long been thought to detect the 
formation of a quark-gluon plasma in ultrarelativistic 
heavy ion collision. Heavy quark radiative energy loss is 
such a probe to study the properties of the quark-gluon 
plasma. In recent years, the investigation of the parton 
energy loss in QCD matter has created considerable 
interest! 1-7 !. The study of the induced gluon radiation 
from heavy quarks in QGP matter is of great importance 
for the understanding of experiment data from high 
energy nucleus- nucleus (AA) collisions at the Relativis- 
tic Heavy Ion Collider (RHIC) and the Large Hadron 
Collider (LHC). Recent measurements' 8-13 ! of high pj_ 
hadron production and its centrality dependence in 
Au-Au collisions at ^/snn = 200 GeV provide the first 
evidence for medium- induced parton energy loss. The 
motivation of this paper is to employ light cone path 
integral (LCPI) approach to solve the massive partons 
radiative energy loss. The analytical solution is obtained 
by the way of light cone path integral. Our analytical 
result is consistent with the former numerical result-the 
parton energy loss tends to be smaller for massive quarks 
than for massless ones. 



In 1953 Landau and Pomeranchuk! 14 ! predicted with 
classical electrodynamics if the formation length of the 
bremsstrahlung becomes comparable to the distance 
over which the multiple scattering becomes important, 
the bremsstrahlung will be suppressed. Migdal' 15 ! 
developed a quantitative theory of this phenomenon. In 
current literature, we call the suppression of radiation 
processes in medium the Landau-Pomeranchuk-Migdal 
(LPM) effect. First results on the LPM effect in QCD 
were given by Gyulassy and wang! 16,17 !, they first 
discussed that the parton jet, produced in inelastic 
collision, propagating in the QCD matter will lose its 
energy due to medium-induce gluon radiation (G-W 
model). They pointed out, comparing with elastic 
energy loss, the contribution of inelastic energy loss 



is more important. G-W model has been extended 
by R.Baier, Y.Dokshitzer, A.H.Mueller, S.Peigne and 
D.Schiff(BDMPS), using equal-time perturbation theory. 
The calculation of BDMPS about the energy loss of 
inelastic scattering indicate that the parton radiative 
energy loss has a square dependence on the path length 
in the medium! 18 !. Based on G-W model, M. Gyulassy, 
P.Levai and I.Vitev(GLV) developed opacity technol- 
ogy to calculate jet energy loss' 19 !. E.K.Wang and 
X.N.Wang first discussed the detailed balance effect of 
jet energy loss in hot QGP medium! 3 !. They pointed 
out the modified energy dependence of the energy 
loss will affect the suppression shape of moderately 
high p t hadrons due to jet quenching in high energy 
heavy ion collision. All above mentioned approaches 
are focused on massless parton energy loss. Based on 
the non-abelian Furry approximation, U.A.Wiedemann 
obtained a medium-induce gluon distribution! 20 ! . It is 
a general result and provides a proof of the equiva- 
lence of the BDMPS and Zakharov formalisms. They 
extended medium-induce gluon distribution to calcu- 
late quenching weig hts [21 l and massive quark energy 
loss' 6 ! and got significative information for jet quenching. 



Due to the mass effect, it is hard to solve the prob- 
lem of massive parton emission. Up to now, most results 
of heavy quark energy loss have been numerical. The 
light-cone path integral approach is a simple way of deal- 
ing with bremsstrahlung of photons and gluons and it 
can give an analytical result. In the path-integral for- 
malism, the radiation cross section is determined by a 
dipole cross section which essentially measures the dif- 
ference between elastic scattering amplitudes of different 
projectile Fock state components as a function of im- 
pact parameter. Baier, Dokshitzer and Schiff (BDMS) 
have shown! 22 ! that the evolution of the rescattering am- 
plitude in the BDMS-formalism is determined by Za- 
kharov's dipole cross section! 20 ! . Accordingly, static De- 
bye screened scattering centers are considered and all 
the scattering centers are supposed to be independent. 
The probability of gluon emission in the LCPI approach 
is expressed through the solution of a two-dimensional 
Schrodinger equation with an imaginary potential. The 
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two-dimensional Hamiltonian reads' 23 ': 

1 , d s9 .n(z)a 3 (p,x) 



H 



2M(x) y dp 



(1) 



where M(x) = Ex(\ — x), x is the gluon fractional 
momentum, n{z) is the number density of the medium 
at the longitudinal coordinate z and 173 is the cross 
section of the interaction of a color singlet qqg system 
with a color center. 

The contribution of the bremsstrahlung mechanism to 
the cross section of gluon production can be written as' 7 ': 



dx 



= Re J dptp*(p, x)a 3 (p,x)ip(p,x,z), 



(2) 



where ip{p 1 x) is the light-cone wave function for the q — > 
qg transition in vacuum and ijj(p,x,z) is the medium- 
modified light-cone wave function for the q — > qg transi- 
tion in medium at the longitudinal coordinate z, p is the 
transverse coordinate and x is the Feynman variable of 
the radiated gluon. The wave functions reads' 7 !: 
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ip(p,x,z) = p(x)( 
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where p(x) = iy / a s /2x [s g (2 — x) + 2s q x]/2M(x), s qg 
denotes parton helicities (s q = 5 , s g = ±1), K is the 
Green function for the two-dimensional Hamiltonian. Kq 
can be written as: 



K (p2,z 2 \pi,z 1 ) 



M{x) 



2-Ki{z 2 - Z\) 



exp[ 



,iM{x){p 2 - pif 



2(z 2 - z x ) 



(5) 

where Lf — 2Ex(l — x)/e 2 ,e 2 = m 2 (l — x) + m 2 x 2 , Lf 
is the gluon formation length. m q is the quark mass and 
m g is the mass of the radiated gluon. The latter plays the 
role of an infrared cut-off removing contributions of the 
long- wave gluon excitations which cannot be treated per- 
turbatively. We assume that the heavy quark is produced 
in the central rapidity region at r\ = and the production 
point is at z = 0, propagating through a hard mechanism 
in a medium of extent L along the z axis. The induced 
gluon bremsstrahlung spectrum can be represented as: 



dp 

dx 



-JZ = / dz <" 



d^f H f (x,z) ^ 
dx 
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where n(z) is the number density of the medium. 

Using the formulas (3)-(5), it is easy to obtain the light- 
cone wave function after a simple calculation: 

p(x)M 2 (x) , 

nP, X ) = ^ (-Px+lSgPy) 
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d£ , if , ,iM(x)p 2 ^ . , 
x /o ^exp(-A)e*p[-L£-]. (8) 

Substituting (7) and (8) to (2), the cross section of gluon 
production can be represented as: 



d<J^f H f(x lZ ) 
dx 



p 2 {x)M 3 (x) 
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with K\ (pe) the modified Bessel function of the second 
kind, 03 (p, x) is the three-body cross section of the imag- 
inary potential ' 24 1 , er 3 (p,x) = Ca/2Cf[& 2 ((1 ~ x )p) + 
(J 2(p) — -^r02(xp)] — C 3 (x)p 2 , where 02 (p) is the dipole 
cross section for scattering of a qq pair on a color cen- 
ter, C z (x) = C 2 (p)A(x) with A(x) = [1 + (1 - x) 2 - 
x 2 /N 2 ]Ca/2C f . In the region p « -j (here p is the 
Debye Screening mass), which dominates the spectrum 
for strong suppression, C 2 (p) takes the form' 7 ! : 



(10) 



where Ca, Cf, Ct arc the color Casimir operators. 
From (10), we can see that C2(p) has a slow logarithmic 
dependence on p. 

In order to derive a quantitative estimate, we take the 
charm quark of mass m q = 1.5 GeV. When pe is small, 
we can expand the Bessel function K\(pe) and keep the 
first two terms. After a complex algebra calculation, one 
can obtain the main contribution of the bremsstrahlung 
to the cross section of gluon production: 
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where G(x) = « s Cf[1 — x + x 2 /2]/x and the Eulcr con- 
stant c ~ 0.5772. The radiative energy loss can be writ- 
ten as: 



AE : 



/ 



dujuj 



dp 
dw 



(12) 



3 



Considering light partons, the mass can be neglected, 
Lf — ► oo, cos -gj w 1, sin w 0. From (11) and (12) we 

obtain the light quark radiative energy loss AE ^ : 



where -J- 



a 2 s TrC F C T A(0)n/2n 2 . At E -> 
loss (13) is equal to GLV's result! 19 !; 

AE GLV = — — In—, 

4 \ g [i 



(13) 

oo the energy 
(14) 



where \i is the Debye Screening mass and 
uj cr <~ max(nC3i 3 /4, L/i 2 /2). Formula (14) reflect 
the logarithm energy dependence of radiative energy 
loss. 

Now let us discuss the radiative energy loss of a heavy 
quark with formula (11). Considering the quark mass 
effect, the gluon formation length Lf is a finite quan- 
tity. In formulas (3), (4), (7), (8) and (9), the exp(-^) 

is no longer equal to 1 as in the case of a light par- 
ton. The mass-dependence of the gluon distribution (6) 
comes from the phase factor exp(— -^), which is anal- 
ogous to the mass dependence of the gluon distribu- 
tion which comes from the phase factor exp[ig(j/; — yi)\, 
q = x 2 m 2 /2u!^\ . In the high energy limit, the formation 
length, Lf, becomes larger than the quark path length in 
the QGP, i.e.Lf » L^l Using (6), (11) and (12), we 
can get: 
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where m g — 0.375 GeV is the mass of the related gluon. 
From (15) we can see that the first term is the radiative 
energy loss of a light quark. Considering the quark mass 
effect, the latter terms are the modifications of the light 
quark. 

For the sake of seeing the suppression of the heavy 
quark energy loss more clearly, we use Fig.l and Fig. 2 to 
show the radiative energy loss of light quarks and heavy 
quarks in a QGP and the dependence on the mass. 
Fig. 1 shows the radiative energy loss of light quarks and 
heavy quarks. In Fig.l we can see that the difference 
between light quarks and heavy quarks becomes smaller 
with increasing energy of the quark. It is obvious that 
when the quark energy is much larger than its mass, the 
heavy quark energy loss is approaching the one of light 
quarks. At high energies, our results are consistent with 
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FIG. 1: The radiative energy loss of light quarks and heavy 
quarks as a function of their energy in a plasma characterized 
by n = 0.5 GeV, A 9 = 1 fm, m g = 0.375 GeV and L = 4 fm. 
The solid line (a) is for light quarks, which were studied by 
Zakharov[7]; the dashed line (6) is for heavy quarks studied 
by us. 
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FIG. 2: The radiative energy loss of light quarks and heavy 
quarks as function of their energy in a plasma characterized by 
At = 0.5 GeV, A g = 1.0 fm, m g = 0.375 GeV and L — 4 fm. 
The solid line (a) is for light quarks, which were studied by 
Zakharov[7]; the dotted line (6) for a quark of m q — 0.2 GeV; 
the dashed line (c) is the energy loss for charm quarks via 
Eq.(15). It is shown AE is dependent on the quark mass, 
the radiative energy loss tends to be smaller for heavy quarks 
than for light quarks. 



Gyulassy's calculations I 25 ' . To illustrate the quark mass 
dependence of AE, we use Fig. 2 to compare the results 
for a heavy quark (charm quark m q = 1.5 GeV) to light 
quarks with mass m q — and m q = 0.2 GeV. With the 
above values of m q , one can see the extent of the AE 
dependence on the quark mass. Both Fig.l and Fig. 2 
show that the formulas (13), (14) and (15) are applicable 
in the high energy limit. At RHIC (^s/vw = 200 GeV) 
and LHC {^snn = 5.5 TeV) energies, all the above 
formulas are operable. 
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In this work, we assume static Debye screened scat- 
tering centers and that all the centers are independent. 
The energy of the parton is supposed to be high enough 
so that the condition Lf >> L is fulfilled. We use the 
light-cone path integral method to deal with the gluon 
emission and thus obtain a simple analytical formula for 
the heavy quark radiative energy loss. Equation (15) can 
be used to estimate the measurable yield of hadrons con- 
taining heavy quarks. To our knowledge, up to now most 
of the calculations of the heavy quark radiative energy 
loss have been carried out by numerical simulations. In 
the high energy limit our results are consistent with Gyu- 
lassy's numerical calculations' 25 !. But one should keep in 



mind that when the quark mass becomes larger, pe is no 
longer a small quantity and the Bessel function Ki(pe) 
cannot be expanded as above. In a following work, we 
will try to extend our result to accommodate all quark 
masses. 
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